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Abstract 

We study the dynamic behavior at high energies of a chain of anharmonic oscil- 
lators coupled at its ends to heat baths at possibly different temperatures. In our 
setup, each oscillator is subject to a homogeneous anharmonic pinning potential 
Vi(.Qi) = ki| 2fc /2fc and harmonic coupling potentials V2(?i — Qi-i) = (n — 1) 2 /2 
between itself and its nearest neighbors. We consider the case k > 1 when the pinning 
potential is stronger then the coupling potential. At high energy, when a large frac- 
tion of the energy is located in the bulk of the chain, breathers appear and block the 
transport of energy through the system, thus slowing its convergence to equilibrium. 

In such a regime, we obtain equations for an effective dynamics by averaging out 
the fast oscillation of the breather. Using this representation and related ideas, we can 
prove a number of results. When the chain is of length three and k > 3/2 we show 
that there exists a unique invariant measure. If k > 2 we further show that the system 
does not relax exponentially fast to this equilibrium by demonstrating that zero is in 
the essential spectrum of the generator of the dynamics. When the chain has five or 
more oscillators and k > 3/2 we show that the generator again has zero in its essential 
spectrum. 

In addition to these rigorous results, a theory is given for the rate of decrease of the 
energy when it is concentrated in one of the oscillators without dissipation. Numerical 
simulations are included which confirm the theory. 

1 Introduction 

One subject that has received considerable attention in recent years is the return to 
equilibrium of systems arising from statistical mechanics. One of the simplest models 
of interest is given by the kinetic Fokker-Planck equation 

d t (f t =£(p t , <po = ip, C= -d%-pd p +pd q -VV(q)d p , p, q € R™ . (1.1) 
This equation describes the evolution of an observable ip : R™ — > R under the Hamil- 

2 

tonian dynamic for the energy H (p, q) = ^- + V(q), perturbed by friction and noise: 



dp = -VV(q) dt~pdt + dw(t) , dq=pdt . (1.2) 
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The relation between (11.2b and (II. lb is given by the fact that the function tptipo, qo) = 
E(fo(p(f), q(t)) satisfies the partial differential equation dl.lt . provided that the pair 
{p(t), q{t)) solves dl.2l ). It can easily be checked by inspection that if V is sufficiently 
smooth and coercive, the measure /.i = exp(— 2iJ(p, q)) dpdq is invariant under this 
dynamic in the sense that if cpt satisfies dl.ll ). then J tptdfi does not depend on t. This 
can be rephrased as saying that if (po, qo) is a random variable with law /i independent 
of the driving noise w, then the law of (p(t), q(t)) is given by fi for all times. 

Underpinning much of the analysis of C on the space L 2 (p) is the guiding principle 
that it related to the corresponding Witten Laplacian 1HN051 

A y = -A, + |W(g)| 2 - (AV)(q) . (1.3) 

In particular, it was conjectured by Helffer and Nier 1HN05I Conjecture 1.2] that C 
has compact resolvent on L 2 (/j,) if and only if Ay has compact resolvent on thenar 
space L 2 (R"). This conjecture has been partially solved in f HN05| (see also QNie06l ) 
in the sense that one can exhibit a large class of potentials for which it holds (loosely, V 
should grow in a sufficiently regular way at infinity). Recently, upper as well as lower 
bounds on the spectral gap of C have been obtained in [HN04 Vil07 1 for potentials V 
that satisfy certain growth conditions at infinity. 

All of the results cited above make heavy use of the following two key facts: 

a. There is an explicit formula for the invariant measure of fll.2t . 

b. The friction term — pd v acts on all (physical) degrees of freedom of the system. 

Both of these facts are violated in the following very popular model for heat con- 
duction. Take a finite collection of N + 1 anharmonic oscillators with nearest-neighbor 
couplings, that is a classical Hamiltonian system with Hamiltonian given by 

N 2 N 

H(p,q) = J2(Y + Vl{qi) ) + ^ V2{qi ' qi - l) ■ (L4) 
i=0 i=l 

Here, the potential V2 is the interaction potential between neighboring oscillators, 
whereas V\ is a pinning potential. This system is then put in contact with two heat 
baths at different temperatures To and T/y. We model the interaction with the heat 
baths by the standard Langevin dynamics, so that the equations of motion of our sys- 
tem are given by 

dpa = -70P0 dt - V{(q ) dt - V 2 \q Q - qi) dt + cr dw Q 
dpi = -V{(qi) dt - V 2 '(<7 4 - q l -i) dt - V 2 \qi - q l+ i) dt (1.5) 
dpN = -JNPN dt - V{(q N ) dt - V 2 \q N - qN-i) dt + a N dw N 
dqj = pj dt . 

Here, we set of = 2%Ti, the index i runs from 1 to N — 1 and the index j runs from 
to N. As described in [BLR00|, the rigorous analysis of this model and in partic- 
ular the derivation of Fourier's law (or the proof of its breakdown) is an outstanding 
mathematical problem of great interest to the applied community. 

If To = Tm = T, then one can check as before that this set of equations has 
a unique invariant measure, which has density exp(— H(p, q)/T) with respect to the 
Lebesgue measure, where H is the Hamiltonian from dl.4b . When the two temperatures 
To and Tn are different however, much less is known. In particular, as we will see 
immediately, even the existence of an invariant probability measure is an open problem 
in some cases as simple as V\(q) = q A and V%(q) = q 2 . 
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This model has been the subject of many studies, both from a numerical and from 
a theoretical perspective. The purely harmonic case has been solved explicitly in 
IIRLL67I . Though no explicit solution is known in the anharmonic case, a wealth of 
numerical experiments exist, see for example IILLP03I and references therein. Since 
we will mainly focus on the theoretical aspects of the model, we refer to HEPR99al 
EPR99b|, which seems to be one of the first rigorous studies of the anharmonic case. 
It was shown in |EH00 EH03] that if V\(q) and V2(q) behave approximately like \q\ ai 
and \q\ a2 respectively at infinity then, provided that a 2 > a\ > 2, there exists a unique 
invariant measure for ( 11.5b . The statement that was proved was actually stronger than 
that, namely it was shown that the generator C of (11.51 has compact resolvent in every 
space oftheformL 2 (exp(-i?(p, q) / T) dp dq) with T > max{T , T N }/2. In HRBT02I . 
it was also shown, using entirely probabilistic rather than functional-analytic tech- 
niques, that the condition a-2 > a\ > 2 is sufficient for the existence and uniqueness 
of an invariant measure for ( 11.5b . Furthermore, the compactness of the corresponding 
semigroup in some weighted L°° space was also proved there. 

This left open the case 02 < ai which is the subject of the present study. To our 
knowledge, no previous rigorous results exist in this case, although some interesting 
theory has been developed recently in HBK07I ILS05I . At first sight, one might think 
that there is no a priori reason for the behavior of ( 11.5b to differ in any essential way 
from the case 02 > a\ where spectral gap results are known. Such wishful thinking 
turns out to be overly optimistic. Even in the simplest possible scenario, that is when 
Vi(.l) = <? 4 an d ^2(9) = <Z 2 , we will show in Theorem l3 . 1 1 I below that the compactness 
property of the resolvent of C is destroyed as soon as N + 1 > 3. Furthermore, when 
N + 1 > 5, it will be shown in Theorem l3.13l that the essential spectrum of C (always 
in a weighted L 2 space of the type considered before) extends all the way to 0. These 
negative results hold even in the case where To = Tn = T, showing that having 
the friction acting on all physical degrees of freedom is a crucial assumption for the 
Helffer-Nier conjecture to hold. 

The reason why the behavior of (11.5b changes so drastically when a.2 < di can 
be understood heuristically by the appearance of breathers (see for example IIMA941 ). 
Breathers are dynamically stable, spatially localized, periodic orbits that arise in the 
noise-free translation invariant (i = —00, . . . , +00) version of (11.5b . Good approxi- 
mations to these orbits persist in ( 11.5b , especially if N is large. It is therefore possible 
to put the system into a state where most of its energy is localized in a few oscillators 
located in the middle of the chain. On the other hand, energy can be dissipated only 
through the terms — ^p appearing in the equations for the first and the last oscillator. 
Therefore, one expects the energy of the system to decay extremely slowly. The appear- 
ance of breathers can be proved in the case where the strength of the nearest-neighbor 
coupling is much weaker than the strength of the pinning potential. At high energies, 
this is precisely the case when 02 < a\ . 

This discussion shows that even the mere existence of an invariant measure is a non- 
trivial problem in this model unless To = Tv, m which case one can check explicitly 
that the usual Boltzmann-Gibbs distribution is invariant. This differs from many other 
systems in that the existence of an invariant measure, not its uniqueness, is difficult. 
In fact, the uniqueness of an invariant measure for a chain of arbitrary length follows 
quickly from the hypoellipticity of the generator and the Hamiltonian structure once 
the existence of an invariant measure is established. 

We are at the moment unable to provide a general proof that shows the existence 
of an invariant measure for a chain of arbitrary length. However, in the case of a chain 
comprising of 3 oscillators, we show in Theorem 15 . 61 below that there exists a unique 
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invariant measure, provided that the coupling potential is harmonic and the pinning 
potential is homogeneous of sufficiently high degree. This result is proven by first 
obtaining an effective dynamics when most of the energy is concentrated in the central 
oscillator. This effective dynamics is then used to construct a Lyapunov function whose 
level sets are compact and hence, by a variation on the classical Kryloff-Bogoliouboff 
method, implies the existence of an invariant measure. 

The remainder of this article is organized as follows. In Section^ we give a formal 
calculation that shows how it is possible to relate the spectral properties of C to the 
scaling properties of the potentials V\ and V 2 . The results given by these formal cal- 
culations are then compared to numerical simulations. We proceed in Section[3]to the 
proof of the negative results concerning the lack of compactness and/or of a spectral 
gap for C. In Section|4]we derive effective equations of motion for the system of three 
oscillators in the regime when a breather is present. Here, we make heavy use of the 
compensator techniques from averaging / homogenization theory [BLP78 1. Finally, we 
show in Section[5]that in the simplest case of a chain of three oscillators with harmonic 
coupling potentials, one can show the existence of an invariant measure for any degree 
of homogeneity greater than 2 for the pinning potential. Perhaps more interesting then 
this last result is the method of proof. We derive as system of effective equations and 
prove their accuracy that when the energy of middle oscillator is large. These effec- 
tive equations give insight into the mechanism of energy dissipation and are used to 
construct a Lyapunov function. 
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2 A formal calculation 

In this section, we first perform a formal calculation that allows us to get a feeling of 
the speed at which energy is extracted from such a system. Consider the simplest non- 
trivial case, that is when N + 1 = 3. In order to keep things simple, we furthermore 
assume that the coupling potential V 2 is quadratic: V 2 (q) = q 2 /2 and that the pinning 
potential is homogeneous of degree 2k: V\{q) = \q\ 2k /(2k) for some real number 
k > 1. The equations of motion for the system of interest are thus given by 

dpo = -7oPo dt - qa\qo\ 2k ~ 2 dt - (q - qi) dt + ^2^T dw 
dpi = -qi \qi\ 2k ~ 2 dt - (2 9l - q - q 2 ) dt 

\2k — 2 I 

dp2 = -I2P2 dt - q 2 \q2 \ dt - (q 2 -q{)dt+ y 2 72 T 2 dw 2 
dqj = pj dt . 

Let us first have a look at the motion of the middle oscillator by itself, i.e. at the solution 
of 

dp = -q\q\ 2k ~ 2 dt , dq = pdt. (2.2) 
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It is easy to see that this equation is invariant under the substitution 

q(t) = E™ q(E a t) , pit) = p(E a t) , (2.3) 

with a = i — Jr. Let us now denote by (p, q) the solution to d2.2l ) with initial condition 
(\/2, 0), so that (p, q) as given by d2.31 l is the (unique up to a phase) solution to ( 12.2I > at 
energy E. Since the variables p and q are assumed to be one-dimensional, the solution 
(p, q) is periodic, say with period r. 

Consider now the equation for the left oscillator, into which we substitute the (ap- 
proximate) solution to the motion of the middle oscillator that we just found: 



dpo = -JoPo dt - q \q \ 2k 2 dt + q dt+ y/2-y T dw + E^ q(E a t) dt 

dq =Podt . (2.4) 

If E is large compared to the typical size of (po, qo) we expect that, up to lower order 
corrections, the solution to this equation behaves like the superposition of the solution 
(Po, Qo) t° (12.4b with the exogenous forcing q = and of a highly oscillatory term of 
the form 

p = E^~ a P(E a t) , P = q. (2.5) 

By symmetry, the same applies to the right oscillator (^2,92)- Applying now Ito's 
formula to the total Hamiltonian H for ( 12. U yields 

^EiT(i) = E( 7o T + l2 T 2 - - I2pl(t)) . (2.6) 

In light of the above discussion, we take (pi, qi) ss (pi + pi,q~i + qi) for i = 0, 2 so 
Pi = Pi + 2piPi + Pi - From its definition, observe that Epf(t) — > Tj as t — > 00 for 
i = 0, 2. Furthermore, when the energy _E of the middle oscillator is large, we expect 
the product po(t)Po(t) to average out to a small quantity when integrated over time 
periods much larger than E~ a due to the highly oscillating, mean zero po(t). Applying 
this line of reasoning to ( 12.6b shows that, in the regime where most of the energy of the 
system is concentrated in the middle oscillator, one expects to have 



^H(t) « -(loVpl + 7 2 Ep|) « -( 7o + 72 ) Kfe i?(i)-- 1 , (2.7) 



d_ 

dt 

where is the variance of the function P introduced in ( 12.5b , that is 



n k = - I P 2 (s) ds-\( f P(s) da) 2 , P(t) = q(t) . (2.8) 



T 

The dependence on fc comes from the fact that q is the position of the free oscillator 
with potential ! -t 2 j-- Actually, one expects this behavior to be correct even in a chain 
with more than just three oscillators. If the chain has N + 1 = 2n + 1 oscillators 
and most of the energy is stored in the middle oscillator, one expects the motion of the 
endpoints to be given by the superposition of a slow motion and a highly oscillatory 
fast motion po with a scaling of the type 



Po 



E^k^ i2n ~ 1)a P(E a t) , 



where P is the (unique) periodic function such that d dt -i n Ji = q and such that the 
integral of P over one period vanishes. (This is because we assume that the nearest- 
neighbor couplings are linear.) This suggests that in the general case of a chain of 
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length N + 1 = 2n + 1, the energy of the system decreases like 

^-H(t) « -(70 + 7Jv)^.n# (i)^ +1 ~ 2n , (2.9) 
at 

for „ the variance of P (of course, Kk,i = ftfc as defined above). If N + 1 = 2n 
is even, the worst-case scenario is obtained by storing most of the energy in one of the 
two middle oscillators. Since their distance to the boundary is the same as the distance 
of the middle oscillator to the boundary in the chain of length 2n + 1, we expect the 
rate of decay of the energy to be similar in both cases. 

Remark 2.1 If the coupling potential is not quadratic but homogeneous of degree 21, 
one can still perform a calculation similar to the one we just did, but one has to be 
more careful. When looking at the influence of qi on pi-\ say, one should take into 
account whether the fast oscillations of qi are of order smaller or larger than 1. If they 
are of order smaller than one, one can linearize the coupling potential. If they are or 
order larger than one, one should multiply them by the scaling exponent arising in the 
coupling. Suppose as before that the chain contains N + 1 = 2n + 1 oscillators and 
that most of it's energy is stored in the middle oscillator (oscillator n). Assume that the 
amplitudes in the fast oscillations of p t and scale like and £^~ Q respectively. 
(Recall that a = | — is the exponent giving the period of the oscillations.) 

One then has (3 n = 1/2. The values of with i < n are given by the following 
recursion formula: 

f {21- l)(/3 m - a) - a if [3 l+1 > a 
ll \ ft+i-2a ifft+i<a 

Using this formula, one can then compute 7 = 2/3 . Note that if £ = 1, one obtains 
(3q = i — 2na which agrees with the value for 7 obtained previously. 

Remark 2.2 One would expect these scaling relations to hold at high energies, even 
if the potentials are not exactly homogeneous. One can then still perform most of 
the analysis presented below by splitting the right-hand side of the equations into a 
homogeneous part and a remainder term and by assuming that the remainder gets small 
(in a suitable sense) at high energies. 

2.1 Numerical simulations 

In this section, we show that there is a surprisingly good agreement, even over ex- 
tremely long time intervals, between numerical simulations of ( 12.11 ) and the predictions 
( 12.71 i and (12.91 l. In order to compare the two, we introduce the function 

H {k > n \p, q) = HP(p, q), (3 = 2n(l - i) . 
Plugging this into ( 12.91 l. we get the prediction 

H {k ' n \t) « H {k - n \0) - (70 + 72 )n(2 - |) Kfc> „ . (2.10) 

A straightforward numerical simulation (essentially integration of the free equation) 
furthermore allows to compute the values of Kfc „ to very high precision. For example, 
we obtain 

K 2 ,i « 0.63546991 , k 3 ,i ~ 0.42363371 . (2.11) 
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We performed numerical simulations for the cases k = 2 and k = 3. Both simulations 
were performed using a modification of the classical Stormer-Verlet method (see for 
example | HLW06 1) to take into account for the friction and the noise. The modification 
was done in such a way that the resulting method is still of order two. 

The simulation for k = 2 was performed at a stepsize h = 10~ 3 , and the simulation 
for k = 3 was performed at a stepsize h = 4 • 10~ 4 . Both simulations used 70 = 72 = 
1.3 and T n = T 2 = 1. 




H{t) 




The prediction obtained from ( 12.101 ) with the values ( 12.1 U is shown as a dashed line 
on these figures, but it fits the numerics so well that it is nearly invisible. We emphasize 
that there were no free parameters in the fit, all constants are predicted by the theory. 
Note that the timescale in the second picture differs by a factor 40 from the timescale 
in the first picture. 

2.2 Comparison to a gradient diffusion 

In this section, we argue that the energy decay rate predicted by ( |2.9t also yields a 
prediction on the qualitative nature of the spectrum of the generator of ( 12. U in the 
weighted space L 2 ip,), where // is the invariant measure. 

The idea is to model the behavior of the energy Hit) by a one-dimensional diffusion 
of the type 

dx = (bix) + a'ix)aix)) dt + aix) dw{t) . (2.12) 
It is well-known that the invariant measure for (12.12b is given by 

Since we expect the invariant measure of ( 12. U to behave roughly like exp(— j3H) dH 
(up to lower-order corrections), we should choose a such that a 2 ix) as \b(x)\ for large 
x. Combining this with ( |2.9t , we obtain the model 

In 

bix) = -x~ l , aix) = x~ l/2 , 7=-r- +l-2n, (2.13) 

k 

which has 2 exp(— 2x) dx as its invariant measure (we restrict ourselves to the half- 
space x > 0). Note that since k > 1 (the pinning potential grows faster than the 
coupling potential), one has always 7 < 1. 

With the choice (12.13b . the generator for (12.12b is then given by 

(Cf)ix) = ^d x (xWJ)ix) - x^d x fix) . 
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Since the operator (Kf)(x) = f(x)exp(—x) is a unitary operator from the weighted 
space L 2 (exp(— 2x) dx) to the fiat L 2 space, the operator C is unitarily equivalent to 
the operator C\ on L 2 given by 

(£i/)(x) = \d x {x~<d x f){x) - ^-(x - 7) . 

At this point, we recall that if ip: R + — > R + is a strictly increasing differentiable 
function with (p(0) = and lim^^oo ip(x) = oo, then the operator 



(U v f)( X ) = f(cp(x))y/^(x) 

is a unitary operator from L 2 (R + ) to itself which furthermore satisfies the identity 
U^ 1 = U v -i. Under conjugation with U v , we see that one has the unitary equivalences 

L :dx - ■ V(x)&V(<p(x)). 



<p'(x) 2(tp'(x)) 2 

Choosing ip{x) = x 2 / (2_7) , we see that C\ is unitarily equivalent to the operator £ 2 
given by 



8 1 I 16 32 

This is a Schrodinger operator with a potential that behaves at infinity like x 2 ' 1 ^ 2 '\ 
It follows that 

• If 1 > 7 > 0, then the operator £ 2 has compact resolvent. 

• If 7 = 0, the operator £ 2 does not have compact resolvent, but it still has a spec- 
tral gap (since one can see that its essential spectrum is the interval [1/2, 00)). 

• If 7 < 0, then belongs to the essential spectrum of C- 2 - 

See for example 1RS78I for a proof. It is then a natural conjecture that the spectrum of 
the generator of ( 11.51 ) on the L 2 space weighted by the invariant measure has the same 
behavior (as a function of the parameter 7 = % + 1 — 2n) as just described. The next 
section is a step towards a proof of this conjecture. 



3 Lack of spectral gap 

The aim of this section is to obtain information on the location of the essential spectrum 
of the generator C for ( 11.5b . This will be accomplished by using ideas from averag- 
ing/homogenization theory to build a set of approximate eigenvectors. Since C is not 
self-adjoint, there are various possible definitions of its essential spectrum (see MEE87I 
or 0GW69I for a survey). We choose to retain the following definition: 

Definition 3.1 For T a closed densely defined operator on a Banach space B, the es- 
sential spectrum a e (T) is defined as the set of all values A g C such that T — X is not 
a semi-Fredholm operator. 
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The set a e (T) is contained in the corresponding sets for all other alternative defini- 
tions of the essential spectrum appearing in the abovementioned works. In this sense, 
the statement "A G er e (T)" used here is the strongest. In particular, it is contained in 
the set 

KEK.(B) 

where JC{B) denotes the ideal of all compact operators on B and er(T) denotes the 
spectrum of an operator T. 

We will make use of the following generalization of Weyl's criterion [RS78 ], which 
gives a useful criterion for identifying the essential spectrum: 

Proposition 3.2 Let T be a closed densely defined operator on a Banach space B. For 
any A G C, A G o~ e (T) if and only if there exist sequences tp n and y>* of elements 
in T>(T) C B and T>(T*) C B* respectively with \\y> n \\ = \\<fn\\ = 1 an d having 
no convergent subsequence such that linin^oo \\Tip n — \(p n \\ = linin^oo ||T*^>* — 
A<|| =0. 

Proof. See I.EE87] Theorem 9.1.3] and IIKat95l Theorem IV.5.1 1], or the original work 
IWol591 . □ 

There are situations in which, even though it is difficult to locate the essential spec- 
trum precisely, one can nevertheless exhibit a sequence <p n as above such that Tip n 
remains bounded. In that case, one has: 

Proposition 3.3 Let T be a closed densely defined operator on a Banach space B with 
non-empty resolvent set. If there exists a sequence ip n of elements in T>(T) C B with 
\\(p n \\ = 1 an d having no convergent subsequence such that lira sup n _ >00 ||T<^„|| < 00, 
then T does not have compact resolvent. 

Proof. This claim follows from the fact that the compactness of the resolvent is equiva- 
lent to the statement that sets of the form {<p | \\ip\\ < 1 & \\T<p\\ < K} are precompact. 

□ 

3.1 Scalings 

Before we turn to the study of the generator, we make some remarks about the regular- 
ity and the behavior of functions that scale in a particular way. 

Denote by Xji. = P8q — Q\Q\ 2k ~ 2 dp the Liouville operator associated to the 
"free" oscillator 

H f {P,Q) = T + ^-. (3.1) 

The constant k is not necessarily an integer, so that in general the function Hf is not 
C°° but only C [2fc] , where [2k] is the integer part of 2k. We introduce the following 
definition: 

Definition 3.4 A function -0 : R 2 \ {0} — ► R is said to scale like if it satisfies the 
relation tf>(XP, A*Q) = X 2a i;(P, Q). 

One has the following elementary result: 
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Lemma 3.5 If ip G C X (R 2 \ {0}) scales like H^, then dpip scales like 2 , dqip 

a-— a+-- — 

scales like Hj 2k , and Xn f ip scales like 2 2k . 

Given some fixed a, there is a one-to-one correspondence between functions cp that 
scale like Hf a and functions on the circle S 1 in the following way. Define a function 
r : S 1 — > R+ by the unique positive solution to 

r 2 {6) cos 2 9 r 2A; ((9)|sin(9| 2fc _ 
2 + 2k ' 

and set 

(Sip)(6) = ip(r(6) cos 9, r(9) sin 6) . 

The function r is bounded away from (it is actually between \[2 and (2fc) 1 / (2fc) ) and, 
by the implicit functions theorem, it is of class C [2fcl . A straightforward but slightly 
tedious calculation shows that, via S, the operator Xu f is conjugated to the differential 
operator 

(X Hf f)(0) = u(6)f(6) , uj(0) = cos 2 9 + r 2k - 2 (9)\ sin 9\ 2k . 
Definition 3.6 A function tp : R 2 — ■> R that scales like Hf a is said to average out to 

** = _ 



oj(9) 

With these preliminaries, the following result is now straightforward: 

Proposition 3.7 Let i[> £ C (R 2 \ {0}) scale like and average out to 0. Then, there 
exists a unique solution if to the equation 

X Hf <p = iP, (3.2) 

q+ j i 

such that ip also averages out to 0. Furthermore, ip scales like 2k 2 and one has 
if G C r '(R 2 \ {0}) with r' = min{[2fc], r + 1}. 

Proof. Let ip be the unique function scaling like Hf a+ " I > ; ~ 5 an( j suc h that 

{ J u(jt) 2n J Q J w(t) 

One can check that one has indeed Xh } <P> = ip- Furthermore, it follows from their 
explicit expressions that both S and S^ 1 map C r functions into C functions as long as 
r < [2k]. □ 

We conclude this section with a small lemma that allows us to compute the L 2 
norm of functions that scale in a certain way. Let (P, Q) G R 2 and (x, y) G R 2 ™ for 
some n > 1. The functions that will be considered in the remainder of this section will 
always be of the form 

<pe(P, Q, x, y) = F(x - g(P, Q), y - HP, QM(P, Q)x(H f ( p > Q)l £ ) - ( 3 - 3 ) 
for some parameter £ > 0. One has 
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Lemma 3.8 Let F e L 2 (R 2n ), let g,h: R 2 — ► R" be measurable, let x- R+ -> R 
be continuous and compactly supported away from 0, and let -0 be a function that is 
continuous away from and scales like for some a € R. Then, the L 2 -norm of tps 
defined as in ( li.il ) satisfies \\feW oc £ Q+ ^+4fc. 

Proof Make the change of variables (P, Q) = (£ip,£^Q) and use the scaling prop- 
erties of ip. □ 

3.2 The case of three oscillators 

Before we tackling the general case of a chain with arbitrary length, let us "cut our 
teeth" on the problem with three oscillators. Since we do not have an explicit expres- 
sion for the invariant measure p, (indeed, at this stage, we do not even know that it ex- 
ists!), we are going to study the generator of dl.51 l in spaces of the type L 2 {e~@ H dp dq). 
As in [EPR99b EHOO ] , it is not expected that the qualitative nature of the spectrum of C 
depends on the choice of j3, as long as (3 < 2 min{/3o, f3 2 } (as usual, we set /3, = 1 /T,). 
Since one expects the true invariant measure to be somehow "in between" the Gibbs 
measures at temperatures To an d T 2 , it is very likely that the qualitative nature of the 
spectrum of C in L 2 (fi) is also the same. 
We write 

= Vl |g2| 2fc <& 
2 2k 2 ' 



H -- 



2 k 



2k 



2 k 



2k 



Remark 3.9 One should not think of Hq (and H2) as the energy one gets by pinning 
qi at 0, but rather as the energy such that the corresponding force is the one that gets 
averaged out over the fast motion of the middle oscillator. 

With this notation, the Liouville operator Xh = d p H d q — d q H d p for the total 
Hamiltonian can be broken up as follows: 

X H = X Ha + X H2 + X Hl + (q + q 2 - 2qi) d Pl + q x (d Po + d P2 ) . 

Recall that the generator of the stochastic dynamics is given by 

C = X H - jop d Po + joT dp - 72P2<9 P2 + -f 2 T 2 d 2 2 . 

The space L 2 (e~^ H dpdq) isometric to L 2 via the operator Kf = e~P H / 2 f. This 
shows that C is conjugate to the operator C = KCK^ 1 on the fiat L 2 space given by 

£ = X H + ^ li T i(i a i ~ a *)Pid Pl + d pi - otiOL*p\ + a-ij 

i=0,2 

= X H + 70T0C + l2 T 2 C 2 0V . 

where we set 

P ,1/3 

a, ■ = — , Ct: = . 

2 1 Ti 2 
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Remark 3.10 Here we see the importance of the condition (3 < 2 min{/3o, P2}' it 
makes sure that the coefficients in front of pf are strictly negative. If this is not the 
case, C is not dissipative anymore and does therefore not generate a Co-semigroup on 
L 2 {e-P H dpdq). 

The main result of this section is: 

Theorem 3.11 Ifk>2, then the operator C does not have compact resolvent for any 
(3 < 2 min{/3o, fo}- Ifk>2, then it has essential spectrum at 0. 

Proof. The aim is to construct a sequence ip n of approximate eigenf unctions such that 
all the ip n are mutually orthogonal, \\ip n \\ = 1 and ||£y>„|| either stays bounded or 
converges to 0. By Propositions [3]2] and [33] this would then immediately imply a lack 
of compactness for the resolvent of C, or even the presence of essential spectrum at 
0. Since the spectral properties of C and of its adjoint C* are the same, we can also 
construct such a sequence of approximate eigenfunctions for C* instead. They can 
then be interpreted as approximate invariant measures for the dynamic ( 11.51 ). Since it 
seems to be a little bit easier to get an intuition about densities of approximate invariant 
measures rather than about approximately invariant observables, this is what we are 
going to do in this section. 
The adjoint of C is given by 

C* = -X H + TblbOO* + 72I2CC)* 
(C ov )* = (a* - a i )p l d Pi + 8%. - aia*p 2 + a* . 

Direct calculation shows that exp(— ctiHi) is an eigenfunction with eigenvalue for 
C l ov and exp(— a* Hi) is an eigenfunction with eigenvalue for (C l ov )* with i £ {0, 2}. 
However, if To 7^ T-2 (and therefore ao ^ 02) it is not possible in general to find 
a closed expression for an eigenfunction with eigenvalue for C*. Note also that 
cti = a* if and only if Ti = 1/(3, which is not surprising since in this case £ l 01J is 
self-adjoint. 

Choose now a function x - R - * [0, 1] which is smooth and compactly supported 
on [1, 2] and set for example £ n = 3", so that the functions x(Hi/£ n ) have disjoint 
support. The formal calculation performed in Section|2]suggests that when the energy 
of the middle oscillator is large, the dynamic of (11.5b keeps that energy approximately 
constant, while the two boundary oscillators equilibrate approximately at Gibbs mea- 
sures at temperatures To and T2 respectively. Our first guess would be therefore to 
build approximate eigenfunctions for C and C* by setting 

tp n = Cne-^-^xiHx/Sn) , <p* = C^o^-^xiH^ £ n ) . 

Here, the constants C n and C* are chosen such that \\ip n \\ — \\<fn\\ — 1- From 

1 1 

Lemma l3~8l one infers that C n oc C* oc £ n 4h 4 . With this guess, we get 

Cip n = C n e a ° n ° 2 2 ^ _ x _ qi ( aopo + a2p2 ) x ^—jj > 

i 

and similarly for C*ip* n . Since p\ f=a En £ n (on the support of ip n ), the first term 
goes to in L 2 . The second term however goes to 00 because of the factor q\, so we 
have to be a little bit more careful in our analysis. 
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The problem is that we have not so far exploited the fact that we also know approx- 
imately what the fast oscillations superimposed over the slow dynamic of the boundary 
oscillators look like, see ( 12.5b . These oscillations can be expressed as a function <I> of 
the middle oscillator, solution to the Poisson equation 

X Hf $(P,Q) = -Q . (3.4) 

By Proposition 13. 71 this equation has a unique solution that averages to along orbits 

j__ i 

of the solutions corresponding to Hf. Furthermore, $ scales like HJ° 2 . In particular, 
we note that $ is bounded when k = 2 and converges to at large energies when 
k > 2. 

Our next guess is therefore to compensate for these fast oscillations by setting pi = 
Pi + (ft) and taking 



■■r<, 



C n exp^-a H (p ,q ) - a 2 H 2 (p 2 , ? 2 ))x(^-) 



Observe at this stage that Hi = i?/(pi, qi), so that we can make use of ( 13.41 when 
computing Xjj^n- One has, for i = 0,2: 

X Hi <p„ = ki| 2fe ~ 2 + qi)®<p n , 

Xh^u = a poqi(pn + a 2 p 2 qnp n , (3.5) 

qi(d Po + d P2 )ip n = —ctopoqiipn ~ a 2 p 2 qi(p n , (3.6) 

I - Pi x' \ 
(So + 92 - 2qi)d Pl (p n = (qo + q 2 - 2gi)( -{a p + a 2 p 2 )dp<& + )ip n , 

^ t"n X ' 

C ov ip n = ai$((ai + a*)pi + ai$)tp n . 

Here, we have omitted the argument (p±, qi) of $ and the argument H\/£ of \ an d x' 
for the sake of simplicity. Note now that ( 13.5b and ( 13.61 cancel each other out exactly. 
It follows from Lemmal3~8lthat 



\(qo + Q2- 2 qi )d Pl ip n \\ < (1 + )(£» 1 + £ n ~ 2 ) < £* 



1 1 



Note that the exponent j + j% appearing in Lemma 13.81 is precisely canceled by the 
normalization constant C n . We also used here the symbol "Pi < ^2 for two expres- 
sions ^i as a shorthand for "there exists a constant C such that *i < C^ 2 " It follows 

j__ 1 

from the above bounds that ||£y> n || < £n 2 ■ 

It is possible to construct approximate eigenfunctions ip* n for £* similarly by setting 



C n evp(-a.oH (p ,qo) - a* 2 H 2 (p 2l q 2 )^jx{^-) 



Note now that the only difference between £ and £* is that one changes the sign of 
Xh and switches ct!j and a*. This shows that the cancellation between ( 13.5b and ( 13.6b 

still takes place when applying £* to so that ||£*y)* || < £ ,f 2 as above. 

If k = 2, it follows that there exists a constant C such that ||£<ys„| + ||£*9?nll < C 
for every 71. If fc > 2, all the exponents appearing the the above expressions are 
negative, so that lim^^oodl^nH + ||£*y>* ||) = 0. Applying Propositions l3.2l and l3.3l 
concludes the proof of the theorem. □ 
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Remark 3.12 It is clear from the proof that the exact same result also holds for a 
chain consisting of 4 oscillators instead of 3. One can construct approximate invariant 
measures in exactly the same way, but one has the additional freedom of choosing to 
take the energy of either of the two middle oscillators to be large. 

3.3 Longer chain 

In this section, we consider a chain of length N + 1 for N > 4. We will show that 
if k > |, then the generator of the dynamic has essential spectrum at 0. To this end, 
define similarly as before 



Hi 



o 



P 2 0+Pj \qo\ 2k + \qi\ 2k (qo-qif + q 2 ! 
2 2k 2 

^ = T+MT+^rJ+.L — 2 — • 

i=3 i=4 



Pi , 1^2 | 2fe 



so that 



= X Ho + X Hr + X Hc + (<7i + q 3 - 2q 2 )d p . 2 + q 2 (d Pl + d P3 ) . 
As in the previous section, we consider the operator C on the flat L 2 space given by 

We have 



Theorem 3.13 If N + 1 > 5 and k > |, f/ien f/te operator C has essential spectrum 
at 0. If k = |, if cioes nof /lave compact resolvent. As previously, these statements are 
independent of the value of (3 < 2 min{/3o, /3at}- 

Proof. Let $ be defined as in the previous subsection. In this section, we do not only 
add a corrector term to pi and p$, but also to qi and q%. We define <E> (2) to be the 
solution to the Poisson equation Xjj f $ <2) ~ $ and we define new variables p and q by 

^ = pi + $(p 2 , 92) , q~i = qi + $ (2) (P2, 92) , 
for 2 = 1,3 and (pi, qi) ~ (pi, qi) otherwise. With this notation, we set as before 
tp n = C n exp(-a H (p, q) - a N H r (p, q))x(H c /£ n ) . 
In order to compute Xn(p n , let us first apply Xh to pi and qf. 

X H qo = Po , 

XhPo = -~q~o\q~o\ 2k ~ 2 + 9i - qo ~ $ (2) , 
X H qi =Pi+ (9o + 92 - 2 qi )d P ¥ 2) , 
XhPi = -9i|9i| 2fe_2 +9o - 2gi 

+ (91 l9i | 2fe " 2 - 91 l9i l 2 ^ 2 ) + (90 + 92 - 2qi)d P <5> + 2$ (2) . 

Hence 



X H H (p, q) = -po$ (2) - (,qMi\ 2k ~ X + 29i - 9o)(9o + 92 - 2 qi )d P ^ 2) (3.7) 
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+ Pi(qiWA 2k ~ 2 - 9iki| 2fe ~ 2 ) + Pi(<Zo + g 2 - 2 qi )d P <S> + 2p 1 & 2) 

= i?i + i?2 + R3 + i?4 + R5 > 

and similarly for XhHAP, q) by symmetry. We have furthermore 

X H H C = p 2 (qi + q 3 - 2q 2 ) , 

so that 

X HX (H c /£ n ) = qi+ J~ 2q2 p2 X '(H c /£ n ) . 
Since C® v (p n = C'o U ip n = 0, we thus have 

Cifin = (—aoX H Ho(p, q) - a N X H H r (p, q) + — ^ttt^tt^ 1 ■ 

V x(Hc/£n) ' 

We bound the terms appearing in this expression in the same way as in the previous 

subsection. Since $ scales like 2 and $ (2) scales like i^ 2fc , Lemma [3~8l shows 
that 

\\RWn\\ <^ _1 , \\RlVnW £ (l+frF)f# _l < £T* , 
\\R3fn\\ <^# _1 , ||H4Vn|| < {l + SF^r 1 <£^~ l , 
ll^nll <£#" X , ||^Vn|| < (l + ^^n 1 . 

-2-_l J__± 

Collecting all these bounds, we obtain ||£^n|| < max{£, 2 ' 1 , f^" 2 }. As before, if 
we set 

ip* n = C n exp(-agi/ (p, q) - a* N H r (p, q))x(H c /£ n ) , 

we obtain the same bounds for ||£*<£* ||. The exponents appearing in all of these bounds 
are strictly negative whenever k > 3/2, thus concluding the proof of Theorem |3.13| 

□ 



4 Effective dynamics 

From now on, we study the case of three oscillators in detail. In this section, we derive 
an effective dynamic for the outer oscillators that is valid in the regime where most of 
the energy is located in the center oscillator. More precisely, we show that there exists 
a change of variable (p,, q{) ^ (p^, q~i) for i = 0, 2 such that the equations of motion 
for (pi, q~i) decouple (to leading order) from the rest of the system, provided that the 
energy of the middle oscillator is large. 

As before, we will use throughout this section the symbol "fi < <3>2 for two ex- 
pressions VfTj as a shorthand for "there exists a constant C such that "Pi < C^-" 
The constant C depends in general on the parameters of the model, but is of course 
independent of the arguments of the 

Theorem 4.1 Assume that k > |. There exist smooth functions & l p and & l q depending 
on (pi,qi,pi,qi) such that if we make the change of variables pi = pi + $^ and 
q_i = Qi + (for i = 0, 2), the equations of motion for (pi 7 q~i) are given by 

dq t = pi dt + Ridt + Si dwi 

_ _ __ _ (4.1) 

dpi = — qi\qi\ 2k 2 dt — qidt — jiPi dt + Rldt + dwi 
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for some adapted processes R l p , R % q , and S* . Furthermore, the error terms R and 
£ satisfy the bounds 

\R l p \ < (E + E 2 )i~ 5 
l^| <1 

/or some 5 > 0. //ere, f/ie energies Ei are given by Ei = 1 + Etfipi, qi). 

Proof. This theorem gives a change of variables where the high-speed oscillations due 
to the presence of a breather located on the middle oscillator (that is the case where 
Eq,E2 -C Ei) have been decoupled from the remaining degrees of freedom, leaving 
an effective "averaged out" dynamic. Recalling the formal calculation performed in 
Section|2l we see that when the energy E is predominantly concentrated in the central 
oscillator, then the amplitudes of the oscillations forp; and qi, i = 0, 2 scale to leading 
order like E*~2 and E*n _1 respectively. This indicates that there are natural break- 
points at k — 2 and k = 3/2. When k > 2, the oscillations of both the pi and the q^ are 
bounded as E increases, so that they can be removed by a change of variables which is 
a bounded perturbation of the identity. 

When k < 2, the amplitude of the oscillations of the pi increases with E, but 
as long as k > 3/2, the amplitude of the does not. This growth will cause extra 
difficulties. If we consider k < 3/2, both amplitudes would grow with E, leading 
to further complications. Since our goal is to outline the ideas without seeking the 
greatest generality, we resist the temptation to analyze all cases and restrict ourselves 
to the case k > 3/2. 

Before we proceed, let us compute the expressions R and £ for a generic choice of 
$p and . Applying Itos formula to pi and q~i, we obtain 

R; = V\ qi + *j) - V'(q t ) + & q + + £% + qi 

4 = £**-**, 52 = ^(1 + ^3*), £* =0-^$*. 

I |2fc 

Here and in the sequel we write V(q) as a shorthand for >%— and oi as a shorthand for 

The case k > 2. The only "bad" term in the equations of motion for (pj, qi) is the 
qi in the right hand side of the equation for pi. The case k > 2 is much easier than the 
case k < 2 since the system is more "rigid" in the later case. One can then simply take 

= and $p = qi), where <i> is the centered solution to the Poisson equation 

X Hf <P = U(P,Q)-Q (4.4) 

and 7Z(P, Q) = Qtf-'(Hf(P, Q)), where ip: R — ► [0, 1] is a smooth function such that 
ip{x) = 1 for \x\ < 1 and ip(x) = for |x| > 2. Making this choice of $p and <I>* in 
( |4~3T > yields 

0, 4 = -*(pi,«i), 

ffi , i?p = -fi$(pi, qi) + d P $(pi,qi)(q + q 2 - 2qi) + TZ(pi,qi) 

Since the function <1> scales like //j fc 2 outside of a compact set, it can be checked 
easily that the bounds ( 14.2b hold, provided that k > 2. Note that in this case, Ei and 



i 




(4.2) 
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Ei are equivalent in the sense that E\ < E{ < Ei since $ is bounded. Therefore all 
occurrences of Ej in the bounds can be replaced by Ej at the cost of multiplicative 
constants. 

The case 3/2 < k < 2. We are now going to assume that k < 2, which the more 
delicate case. Note that the second and third terms in R l p above satisfy the bounds d4.2b 
(with Ej instead of Ej, but this problem will be dealt with later) provided that QdpQ 
scales like H e for some 9 < 0. This is the case when k > §, which is one of the 
reasons why we restrict ourselves to this case. Therefore, only the terms involving $ 
scale worse then the desired bounds on the error terms and need to be eliminated. This 
motivates the introduction of the solution $ (2) to the Poisson equation 

X ff/ $ (2) = $. (4.5) 

Note that $ (2) scales like H^^ 1 outside of a compact set, so that it is bounded if 
k > |. It would be tempting at this point to simply subtract 7i ( f )(2) (pi, Qi) to pi and 
add $ (2) (pi,gi)to ft . 

This would however introduce correction terms that grow faster than the bounds in 
(I4.21 l. The trick is to realize that these correction terms are multiplied with terms that 
go to as the energy of the middle oscillator becomes large. We therefore multiply <I> (2) 
with a cutoff function that makes sure that this second round of correction is applied 
only when the energy of the middle oscillator is large compared to the energy of the 
boundary oscillators. 

Let ip: R+ — > [0, 1] be a smooth increasing function such that ip(x) = for x < 1 
and ip(x) = 1 for x > 2. Let furthermore a be an exponent to be determined later and 
set 

% = Ql ) - 7i $ (2, (pi, qx) , <E2o 
& q = <p(Ei/E? )$ (2) (pi, qi) , &) 

where we defined £7, = 1 + Hf{pi, qi). In the sequel, we are going to use the shorthand 
(pf = (p(Ei /Ef) and we will omit the arguments of $, $ (2) and cpf in order to simplify 
notations. Before we turn to the verification of the bounds ( 14.21 . we remark that since 
k < 2, Ei and Ei are not equivalent for i = 0,2. Since k > 3/2, E\ and E\ are 
however equivalent in the sense that E\ < E\ < Ei. Since we wish to bound the 
remainder in terms of powers of Ei and not Ei, we are now going to show how these 
quantities are related. From the definitions of pi and qi, we have for i = 0, 2 the 
estimate 

\Ei - Ei\ < 1 + < 1 + epf + < 1 + sEi + -E?' 1 . 

By choosing e sufficiently small and moving the term sEi to the left hand side, we thus 
obtain the two bounds 

Ei < Ei + E^ 1 , Ei < Ei + E^ 1 . (4.7) 

Writing g = E\/Ef, and applying Ito's formula to g, we obtain 

2 

dg = Er a Pl (q + q 2 - 2 Ql ) dt - aE^- 1 ^ - 7?? , 2 + ?1 - q t ) dt 
+ ^^ E^r^olvl dt - aE^r—^iPi dw t 
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= hdt + I 2 dt + I a dt + dM(t) . 

M is a martingale since, on any finite time interval, it is easy to get control over the 
expected value of any power of the total energy. (Even a bound which grows expo- 
nentially with the length of the time interval is sufficient.) Bounding each of the three 
terms Ij separately we obtain: 

\h\ < E; a+ ^E? + E^Ef^El + Ef a E? + ^ , 
\I 2 \ < EiE'" + El + * E^- 1 , 
|7 3 | < BkE^- 1 • 

Applying Ito's formula to ipf, we get 

C<p? = ip'{E x /Ef)Cg + a 2 ^"(E 1 /E^ElEr^ P j = J 4 +I 5 . 

Note now that ip' and ip" are zero outside the interval [1,2]. Therefore, we have E\ < 
Ef < E\ on the support of these functions. This line of reasoning yields for 1/fc < 
a < 2k the bounds: 

\h\ < Ef (i ~ a) + E^Ef^ + E?^ + 1 + Ef- 1 + E^ < , 
\h\<E^<l. 

Collecting these estimates and taking into account the support of ip' and ip" produces 

\£^\<(l + Ef^)l 2E? > El > Er . 

Recall that, from the scalings given in Lemma [331 of the solutions to the Poisson 
equation ( 13.2b and the definitions ( 14.41 ) and ( 14.5b of $ and <i> {2) , we have that |$| < 

Ef~K \d P $\ < ~\ |$ (2) | < E^' 1 , \d P ¥ 2) \ < Ep~^. We will also use the 
fact that, from the definition of <pf combined with ( 14.71 ), one has both bounds 

Ef < E x , Ef < E x , (4.8) 

on the support of ipf. Furthermore, one has the bound Ei < Ef on the support of 
1 — ipf which, by ( |4.7t , implies that one also has E x < Ef on the support of 1 — ipf. 
Using the definitions of and given in ( |4.6t , equation d4.3l ) yields 

+ ipf£& 2> + (&pf)¥ 2) + 7l $ (2) • 

We now make use of the definitions of C and <1> (2) to obtain 

£$(2> = $ + ( qo + q . 2 _ 2q 1 )d P ^ {2) . 

This allows us to obtain the following bounds for R l q : 

\R\\ = \(cpf - 1)$ + ipf(q + q 2 - 2gi) d P ¥ 2) + (£^)$ (2) + 7 ,<i> (2) | , 
< (1 - <p?)Ef~* + ipfE~ ¥l ~^{Ep + Ep + Ejtf) 

+ E^E^Ue^e^ec. + 1 
= J 6 + I r + h + 1 • 
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Our aim is to bound the terms I§, It, 1% in terms of the E/s instead of the Ej's. To do 
so, we now fix 

3 

a= 2- 

Since k £ (3/2, 2), a £ (1/fc, 2fc) which is the constraint that we had to impose earlier. 
Since, as mentioned above, E\ < Ef on the support of 1 — cpf , one has 

\h\ < Et f ■ 

One can check that the choice a = 3/2 implies that there exists 5 > so that f — % < 
jj: — 5 (actually, one can set S — 1/12 for the range of values of k considered here). 

_ J. c 

This shows that |/ 6 | < Ef h as required. 
To bound Ij, we make use of the fact that 

Ep + < Ep + EY~^ + Ep z < Ep + Ep + Ep t , (4.9) 
so that, by virtue of ( 14.81 ). 

\M < <p?E; i(1 -*\Ep + Ep t ) + ^e\-^ 

< E7 s (Ep +Ep l ) + 1 . 

We finally turn to 7 8 . We first remark that, by ( 14.71 ) and the fact that a(| — 1) < 1, 
one has Ef < E\ < Ef on the support of the indicator function 1a where A — 
{2Ef >E 1 >E?}. We thus obtain 

\h\ < Ep^Ef- a \ A < (Ep t + Ep^-^Ep^U 

-J- --2.c2._ll -J_ - i - a 1 + k - k2 

<(E™ i + E^ < E^E-o + E l . 

One can check that, for a = 3/2 and the range of /c's of interest, the exponent of 

-_ x -j_ _j_ 

the last term is strictly smaller than l/2fc. Therefore, |I 8 | < E i {Ef + as 
required, choosing S smaller if necessary. Collecting all of these bounds shows that R l q 
does indeed satisfy the bound advertised in ( 14.2b . 
Turning to R l p , ( 14.31 yields 

141 = 1^(90 + (*> + Q2- 2<?i) (_ P _■ - 7 ^p$ (2) ) - (7? - <^)* (2) 
+ V'{ qi + tf& 2) )-V'{ qi )\, 

< 1 + (fj* + Ep + Ep)E^ 1 + 1 + \V'(q t + ^$ (2) ) - V'( qi )\ . 

Here, we made use of ( 14.91 ) and of (14.7b . Since, for fc > 3/2, $ (2) is a bounded function, 
one has the further bound 

\V% + tf& 2) ) - V'(qi)\ = \V'(qi) - V'(qi - ^$ (2) )| 

< (1 + |*|) 2fe " 2 |$ <2) | < E]'~ k < Ef~ S , 

for 8 sufficiently small. Collecting these bounds and using the fact that k > 3/2, we 
obtain 

\K\ < E ^ + E * + E t S Z Ef~ 5 + Et! , (4.10) 
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which is indeed of the form (14.21 ). 

Lastly, the S -terms can be bounded by 

|E*| = \aa i <p'(E 1 /E?)&»E 1 Er«-i pi \ < > 

where we made use of the fact that, like in the bound of Ig,, E\, Ef and Ef are 
equivalent on the support of ip'(Ei/ Ef). This concludes the proof of Theorem l4.ll 

□ 

In the case k > 2, it will be useful in the sequel to have a better approximation of 
the dynamics that yields smaller error terms in the regime where most of the energy is 
located in the middle oscillator: 



Theorem 4.2 Assume that k > 2. There exist smooth functions Q p and depending 
on (pi,qi,pi,qi) such that if we make the change of variables pi = pi + & p and 
Qi = Qi + (far i — 0, 2), the equations of motion for (pi, qf) are given by 

dcji = fii dt + R l „ dt 

Z (4.11) 

dpi = — q~i\q~i\ 2k ~ 2 dt — q~i dt — ^ipi dt + K l p dt + er,; dwi 

for some adapted processes R p and R l q . Furthermore, the error terms R satisfy the 
bounds 

\K\ < (Eo + ErfH^- 1 \R q \ < (E + E 2 )H^- 1 . (4.12) 

Here, the energies E^ are given as before by Ei = 1 + Hf(pi, qi) and H is the total 
Hamiltonian of our system. 



Proof. Following the proof of Theorem l4.ll we set as in (14. 6t 

% = $(Pi, Qi) - li& 2) (Pi, <?i) . 

$; = $ (2) (pi,gi)- 

This yields for R l p and R q the expressions 

R\ = (qo + q2 - 2 qi )d P ¥ 2) + 7l $ (2) , 

R p = TZ( qi ) + (q +q 2 - 2qO (dp^> - ^d P & 2) ) - ( 7 2 - 1)$ (2) 
+ V'( qi + $ (2) ) - V%) . 

The desired bounds now follow from Lemma POl below. together with the fact that both 
$ (2) and Qd P <$> scale like H^' 1 . □ 

Lemma 4.3 For every a > 0, the bound 

x~ a <(x + y y a max{2 Q , y a } 

holds for every x, y > 1. 

Proof. lfx>y, then x~ a < 2 a (2x)~ a < 2 a (x + y)~ a . If on the other hand x < y, 
then x~ a < 1 < y a (x + y)~ a . □ 
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5 Existence of an invariant measure 
5.1 General strategy 

To prove the existence of an invariant measure, our aim is to construct a type of Lya- 
punov function V{p, q) such that V(p, q) — > oo as \(p, q)\ — ► oo and CV(p, q) — » — oo as 
\(Pi Q)\ ~ ¥ oo. Given such a function, the existence of an invariant measure will follow 
from the following proposition which is a variant of the classical Kryloff-Bogoliouboff 
construction HKB37llHas80l p. 52]: 

Proposition 5.1 Consider an SDE on R" with smooth coefficients and denote its gen- 
erator by C. Assume that the SDE has global solutions and generates a Feller semi- 
group. If there exists a smooth function V: R" — > [0,oo) such that the level sets 
{x : CV(x) > C} are compact for every C, then the SDE possesses an invariant 
probability measure p. Furthermore, the function CV is integrable with respect to p 
and J CV(x) p(dx) = 0. 

We will construct the function V in steps by analyzing the dynamic in the limit 
of various energies being large and then draw inspiration from the structure of these 
limiting regimes to construct V. Operationally, we will make an initial guess for V and 
then augment it by a series of correction terms. 

Proof of Proposition ^. 1\ Denote by x t the solution to the SDE starting at some (de- 
terministic) initial condition xq. Applying Ito's formula to V(xt), we get 

dV(x t ) = (CV)(x t ) ds + dM(t) , 

for some continuous local martingale M. Therefore, there exists an increasing se- 
quence Tjy of stopping times converging to +oo such that M(t A rjv) are martingales. 

EV(x tATjv ) - V(rro) - E / (CV)(x s ) ds = . 
Jo 

Since V is positive this shows that, for every K > 0, 

E / (K — (£V)(a? 8 )) ds < V(x ) + Kt . 
Jo 

Taking K large enough so that K — CV > 0, we can apply the monotone convergence 
theorem to take the limit TV — > oo and obtain 

-E f {CV)(x s )ds <V{x ) . 
Jo 

Now, by assumption the sets Ar = {x : —CV(x) < R} are compact for all R. In 
particular, this implies that there exist a K > so that — CV(x) + K > for all x. 
Now observe that for R > —K 

- / P(a; s g A R )ds = - I P( - CV(x s ) + K > R + K)ds 
t Jo t Jo 

< 1 r* K-KCV(x„) < K + V(x ) 

- t J R + K S ~ K + R 
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Therefore, the sequence of measures /j, t defined for measurable sets A by 



1 



Ht(A) = - / P X0 (x s eA)ds 
1 Jo 

is tight. Hence the Kryloff-Bogoliouboff construction |KB37l lHas80l p. 52] guarantees 
the existence of an invariant measure. The last statement then follows from Lebesgue's 
dominated convergence theorem. □ 

Remark 5.2 We will actually be able to construct a positive smooth function V with 
compact level sets that has the property that 

CV < C\ - C 2 V a , 

for some positive constants Cj and some (typically quite small) a <G (0, 1]. In this 
case, it is known MDFG06I that one does not only have the existence of an invariant 
measure, but the transition probabilities converge towards it at rate 0(t~ Q / (1 ~ Q) ) in the 
total variation distance. We believe that this convergence actually takes place at a much 
faster rate, but such a statement is beyond our reach at the moment. See also lVer97l 
for related results on subexponential mixing for SDEs. 

5.2 Construction of the Lyapunov function 

Recall the change of variables (p, q) i— » (p, q) from Theorem l4. 1 I that leads to an effec- 
tive decoupled dynamic for the outside oscillators. In order to construct the Lyapunov 
function V, we proceed in two steps: 

1. We gain good control over the dissipation of the energy stored in the outside 
oscillators. This will be the content of Proposition l5.3l 

2. We use this in order to get control over the dissipation of the total energy of the 
system in Theorem l5.6l 

Proposition 5.3 There is a function Uq equivalent to Hf(pQ, qo) + Hf(p2, q 2 ) and such 
that, for every m > 0, there exist constants C m > and c rn > such that CUJf 1 < 
C — r 7/ m 

Proof. Inspired by ( 14. U . we define an effective Hamiltonian Hq by 

(Note that Hq is equivalent to Hf.) We set 

Uo = HoiPo, Qo) + Ho(P2, 92) + 7(Po<Zo + 

for some constant 7 to be determined later. If 7 is sufficiently small and since k > 1, 
this function is indeed equivalent to Eq + E 2 — Hf(po, qo) + Hf(p2,q~2)- Applying 
Ito's formula to it, we get from ( 14.11 ) that 



i=0,2 



dUo=J2 ( ( T " 7i)Pi ~ 7(lft| 2fc + l%| 2 ) - lliMi + (Pi + 7?i)i$ 

m 2 "" 2 ^ + iP i )B i q + (k - i)!^ 2 *- 2 !^! 2 + ||s;i 2 + 7 s;s;) dt 
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+ J2 ((P* + ^) S p + (l«l 2 * -2 * + « + 7Pi)Sj) ■ 

i=0,2 

Fixing 7=5 min{l, 70, 72}, we obtain for some constant C 

dU < -2^U dt + C (\PiQi\ + \Pi+7qi\\^ P \ (5.1) 

t=0,2 

+ (Iftl 3 *- 1 + 1*1)141 + (i + Iftl 2 *" 2 )^*! 2 + \Kf) <tt + dM(t) . 

Here, M is a continuous Martingale with quadratic variation bounded by 



Using the notation and the bounds of Theorem |4.1| we have 

\m\<E? + ", \Pi + lQi\\K\^E]~ S 



\U 2k - X + \Pi\)\R\\ <> E\- S + Elzf , (1 + \U 2k ~ 2 )\K\ S E : 



2 fc 



Since all the powers appearing on the right hand sides of these bounds are strictly less 
than 1, we have shown that 

CU < C - 7W0 , 

for some (different) constant C. 

In order to bound the quadratic variation of M, note that one has 

Op.I 2 + |gi| 2 )i4i 2 < Si , (1 + m ik - 2 + |*| 2 )|s;i 2 < e]- l * . 

In particular, one has for some positive constant C 

(M)(t) <C [ (1+U (s))ds . 
Jo 

The claim now follows by applying Ito's formula to U™. □ 

Remark 5.4 Observe that one could also apply Ito's formula to exp(9Uo) for suffi- 
ciently small 9 and obtain C exp(9Uo) < C — a exp{9Uo). 

Remark 5.5 If k > 2, then Uq is also equivalent to Hf(p , q ) + Hf(p 2 , 92)- This is 
not the case when k < 2. 

We are now ready to prove the main theorem of this section. Recall that H is the 
total Hamiltonian defined in ( 11.4b . 



Theorem 5.6 Consider the equations of motion ( 12.71 ) with k > 3/2. Then there exists 
a function V and constants c, C > such that V > cH a — C and such that CV < 
C — cH a for some exponents a and a'. 
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Proof. The idea is to work "modulo powers of Uq." Assume that we can find a function 
U\ such that 

Ui > cH a - CU^ , CUi < -cH a ' + CU^' , (5.2) 

for some positive exponents a, a' and some (possibly very large) exponents N, N'. 
We claim that it then suffices to take V = U x + U^ +N ' +1 . Note that 

V>cH a - CU* + U^ +N ' +1 > cH a - C , 

for some constant C , and so V grows at infinity (has compact level sets). It then 
follows from Proposition 15.31 and from ( 15.21 ) that 

CV < -cH a ' + CUg' + C N+N , +1 - c N+N , +1 U* +N ' +l < C" - cH a ' , 

for some constant C", which is the desired bound. It therefore remains to construct a 
functional satisfying ( 15.21 ). 

The starting point for the construction of U\ is to apply Ito's formula to H n , Since 

dH = J2i=o,2(-7iPi dt + |crf dt + p t Ui dwi), we have 

2 

dH n = nH n - 2 (#(y - Hpf) dt + \{n - l)afpf dt + Ha lPi dw^j . (5.3) 

1=0,2 

At this stage, we have to distinguish two cases, as in the proof of Theorem l4.ll 

• If 3/2 < k < 2, then the motion of the interior oscillator induces large (i.e. 
going to infinity with the energy of the middle oscillator) fluctuations in the 
values of po and p2- Therefore, there will always be energy dissipation. 

• If k > 2, then the motion of the interior oscillator induces bounded (or small) 
fluctuations in the values of p n and p 2 . In this case, one can subtract a compen- 
sator so that only fluctuations remain. There will be a few error terms that seem 
to be larger than the dominant one, but they can hopefully just be eliminated 
order by order. 

The case 3/2 < k < 2: In contrast to the arguments in the proof of Theorem l4.ll this 
is the "easy" case for this part of the proof. It follows from the proof of Theorem 14. II 
that in this case pf is equivalent to pf with pi = pi — $(pi,<7i). Since 2jipi$ < 
■T^ji® 2 + 2jipf, this allows us to obtain for CH n the bound 

CH n < nH n -\C n + j p 2 + l2P 2 - ±(7o + 7 2 )$ 2 ) , 
for some positive constant C n . This in turn implies that 

£H n < H n-l {(Ji + C2Uo _ | (7o + 72)$2) > (5 4) 

for some constants C[. Observe that one can use a similar calculation to obtain the 
bound 

\CH n \ < H^id + C 2 U + C 3 $ 2 ) , (5.5) 

for some possibly different constants Cj. 

Note that $ 2 scales like HI 1 which goes to infinity at high energies. We can thus 
find a positive constant Kk (which was introduced in (I2.81 l) and a function 1Z' such that 
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<I> 2 — KkHJ 1 + IZ'(Hf) averages out to in the sense of Proposition [321 We define 
to be the solution to the Poisson equation 

X Hf V = $ 2 - KkHf' 1 + K'(H f ) . (5.6) 
This function is smooth since TV is precisely such that $ = in a neighborhood of 

_5 3 

the origin. The function $ then scales like Tip 2 . In particular, and Qdp^ are 
bounded. 

With all these preliminaries done, we define 

Ui =H" + f(7o + 72) J ff' 1 " 1 * • 
Applying Ito's formula to U\ and using ( 15.41 ), ( 15.51 , and ( 15.61 ). we obtain 

CUi < CH n - 1 (l+U ) + C\^\H n - 2 (l+Uo + ^ 2 ) 

+ f (70 + 72)ff"- 1 (-«^i l_1 + (90 + Q2- 2qi)d P V) . 
Since |>I>|<I> 2 < H and |(g + 92 - 2qi)dp^\ < U , we have 

CUx < -^(7o + ^W^KuEf- 1 + CH n ~ l U , (5.7) 

for some constant C. 

We can check that U\ satisfies the first bound in ( 15.2b . since \I> grows slower than TT 
and Uq is bounded from below by a positive constant. However, it is not clear a priori 
from this expression that it IA\ satisfies the second bound in ( 15.21 ). 

One can check that 

U a > c(E + E 2 ) - C<I> 2 , 

— — l 

for some positive constants c and C. Furthermore, <I> 2 scales like Ti^ which is 
strictly smaller than Hf, so that 

Uq + E x > cH . 
In particular, there exists a constant C such that 

Hi- 1 < C{Ef^ . 
Inserting this into ( 15.7b . we obtain 

CUi < -cH n+ i- 2 + CTT n - x U a , 

for some constants c and C. Since f — 1 > 0, this shows the existence of a (sufficiently 
large) power N such that 

CUi < -cH n+ i~ 2 + CUg , 

which is the desired bound. 

The case k > 2: The reasoning above only worked for k < 2 since in this case one has 
n + | — 2 > n — 1. When k > 2 the situation is slightly more delicate. Here, as before, 
the energy dissipation mechanism comes from fluctuations of p\ and p\ around their 
mean. However, the amplitude of these oscillations now decreases as the energy stored 
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in the middle oscillator increases, so that we have to be more careful. In particular, p 2 
can no longer be treated as a perturbation with respect to $ 2 , as we did in the previous 
case. 

The expression in ( 15.3b suggests that, in order to extract these small fluctuations, 
we should compensate H n by subtracting 

5 = nH'^lHoipo, g ) + H (p 2 , &)] , 

where the "bar" variables are as in Theorem 14.21 Note that here, and anywhere from 
this point on, we use the variables (pi, q~i) from Theorem l4.2l and not from Theorem l4.ll 
This is because we require a sufficiently good effective dynamics so that the error terms 
are small with respect to $ 2 . The bounds obtained in Theorem l4.1l are not sufficiently 
small for that. Setting as before Ei as a shorthand for Hf{pi, q~i), one has 

CE = nH n - 2 (#(y - 7iP?) +(n~ D^(f - 7iP?) 

i=0,2 

+ i(n - l)(n - 2)E l H~ 1 a 2 p 2 + (n - l)<r?p#i 
+ i?fei?; + ^(|g,| 2fc - 2 + l)i?^)) 

The important term in this expression is nH n ~ x { — 9 g2 — 70P 2 , — 72P2)' a ^ me other 
terms will be treated as perturbations. Setting U\ = H n — E, we obtain 

£U° = -nH n - x ~ Pi) + £ ■ 

j=0,2 

with 

£ = nH n - 2 (-(« - y - 7iP? ) - H{piK + Vi(Wi\ 2k ~ 2 + ^4) 

i=0,2 

- |(n - l)a 2 ((7 1 - 2)E l H- 1 p 2 - 2pm + p. 2 )) . 

We can check that one has \£\ < i7' l ~ 2+ A^2 so that, for every S > 0, there exist 
constants C and N such that the error term £ is bounded by H n+ ik~ 2+d + CUq . 
At this stage, we use the explicit expression for pi to rewrite this bound as 

C1A\ < -nH' 1 - 1 Y 7»(* 2 + 7 2 2 (* (2) ) 2 - 2 7i $$ (2) + 2p t ( lt ¥ 2) -$))+£. 

i=0,2 

This puts us now in a situation similar to d5.4t , with the difference that, up to powers of 
Uq, the error term £ is of order H n ~ 2+ 2k + s instead of being of order H n ~ x . Since $ 2 

is larger than Hp , this is the feature that will allow us to obtain the required bound. 

Since $ (2) scales like Hj k , we obtain in the same way as in the proof of Theo- 
remgjthat |$ (2) | < H A" 1 ^ 1- *. This shows that 

CU\ < -nH' 1 - 1 Y 7($ 2 - 2ft$) + £ 2 , (5.8) 

i=0,2 

where the error term £2 satisfies the bound ^ H n+ &~ 2+s + CUq as befor^H 

1 The interested reader can check that expanding the square around pi instead of pi would lead to trou- 
blesome terms. 
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Since these terms oscillate very rapidly, we would like to replace them by their 

averaged effect over one period of the middle oscillator. To leading order, the terms 

2 — 1 

Pi$ and $ 2 will average out to and k^E* respectively. The latter contribution will 
turn out to be the dominant term leading to an overall dissipation of energy. Defining 
and "J*' 2 ' as in ( 15.6b and d4.5l ) respectively, we finally set 

U x = U{ + nH n ~ 1 ^2 - 2 K $(2) ) ■ ( 5 - 9 ) 

i=0,2 

Our investigation of CU\ begins with 

CH n -H = -H n - x ( Kk Ef- 1 - $ 2 ) + W 1 - 1 (n'iEx) + (q + q 2 - 2 gi )d P ^ 

2 

+ (n - l)H n ~H ( H (-Y - 7<P?) + ¥ n - V°lvi) 

i=0,2 

= -U n - x {K k E\- X - $ 2 ) + h + ho ■ 

Since the function 1Z' has compact support, we have for example \TL'{Ei)\ < E^ 1 . 
Using the trick from Lemma |4~3l this yields \R!(Ei)\ < H^Uq. This allows us to 
obtain the bound 

\h\ < H n ~ 2 U + H n - 1+ ^Ep~ 2 < H n - 2 U + H n - 3+ ^U^* . 
Since \I> is bounded, one furthermore has 

\ho\ < H n ~ 2 U + H n - 3 U . 
Turning to the second term in OA\ , we have 

CH n - l p^ {2) = H^Pi® + H n - 1 pi(q + q 2 - 2 qi )d P & 2) 

2 

+ (n - l)H n - 3 p^ 2) Y, (#(y ~ ljP ^ + 5 (n _ 2)f7 ^) 

3=0,2 

+ H n - 1 ^ 2 \-q i \q i \ 2k ~ 2 - ft + 4 - liPi) +(n- l)H n - 2 a 2 Pl & 2) 

= H^Pi® + In + I 12 + Il3 + hi ■ 

Using the bound on R l p from J4. 12b . the error terms appearing in this expression can be 
bounded by 

I Jul ^ H n ~ 1+ ^U^E^^ < H n -^+'iul~^ , 
\h2 1 < H n ~ 2 U^ + H n ~ 3 U^ , 

\hs\ < H^Ep^iUv+UlH^- 1 ) < H n - 2+ ^U^* , 
\I U \ <H n - 2 ul . 

By Young's inequality, it can be checked that, for every S > there exist constants C 
and N such that \Ij \ < H n - 2+ ^ +s + CUg for j = 9, . . . , 14. Inserting these bounds 
into (l5T8T > and (IBT^I l, this shows that 

OAi < -nK k ( 7o + 72 )/J"- 1 /5*" 1 + H n - 2+ ^ +s + CU* . 
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Since § — 1 < and since n — 2 + ^ >n — 2 + for sufficiently small S, we 

finally obtain 

for a possibly different constant C. This is precisely the bound ( 15.2b which was the 
missing piece to complete the proof of Theorem l5.6l the principal result of this article. 

□ 
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